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I Abstract 

We study single top production at the LHC in a SUSY-QCD model with a heavy 
Dirac gluino. The presence of a heavy Dirac gluino allows for notable top-up flavour 
changing neutral currents. In this scenario, we find that the process ug — > tg gives 
the largest contribution to single top production via FCNCs at the LHC. The key 
features of this signal are that the top quark is produced very forward and that it 
is asymmetric to its anti-top counterpart, as the latter lacks a valence quark. 



1 Introduction 

Due to its large mass, the top quark's presence as an initial partonic state at hadron 
colliders is negligible. The single production of top quarks must, therefore, proceed via 
flavour changing interactions. In the Standard Model (SM) the tree level couplings of the 
W bosons to quarks generate all such interactions. At the loop level, flavour changing 
neutral current (FCNC) interactions are possible. Single top production via FCNCs in the 
SM is, however, strongly suppressed by the Glashow-Iliopoulos-Maiani (GIM) mechanism 
due to the small mass differences of the down type quarks occurring in the loop. 

In supersymmetric extensions of the SM, flavour mixing in the squark sector allows for 
additional top quark FCNCs. Unlike the down type quarks, the mediating squarks may 
have suitably different masses to avoid a similar supression from a GIM-like mechanism. 
In the Minimal Supersymmetric Standard Model (MSSM), however, squark flavour mixing 
is strongly constrained due to the excellent agreement between the SM and experiment [H 
[21 [3j IH [5] . The flavour mixing of the first generation of squarks with those of the second 
and third are constrained by the K°, B% and D° neutral meson mixing experiments. 
Likewise, the flavour mixing of the second and third generations of squarks is constrained 
by B° s mixing as well as experimental results from the FCNC process h — > s'j. 

Phenomenological studies of top quark FCNC processes in the MSSM have predom- 
inantly focused on the less-constrained stop-scharm mixing [61 [7J [91 [101 El 112] ■ The 
dominant contributions are found to come from the SUSY-QCD sector, namely, from 
gluino, as opposed to neutralino, exchange. For a review of top quark FCNCs in new 
physics and a more detailed list of references see, for example, Ref. [T3] . 

Recently, it has been shown that sizable squark flavour mixing is possible in models 
with heavy Dirac gauginos [TH 15J. Specifically, it was found that heavy Dirac gauginos 
suppress the SUSY contributions to neutral meson mixing, thereby relaxing the associated 
experimental constraints on squark mixing. The Dirac gauginos were considered in the 
context of the Minimal R-symmetric Supersymmetric Standard Model (MRSSM): a SUSY 
extension of the SM in which the global U(1)r symmetry of the N=l superalgebra remains 
unbroken. In such models gauginos are required to be Dirac fermions because their charge 
under the continuous R symmetry forbids a Majorana mass term. One means of promoting 
gauginos to Dirac fermions involves the addition of adjoint chiral superfields together with 
a supersoft breaking mechanism [16] . Alternatively, Dirac gauginos appear automatically 
in the hypermultiplets of N=2 SUSY. 

In this paper we study the effects of the extra squark mixing freedom that heavy 
Dirac gluinos allow on single top production. In this context the squark mixing of the 
third generation with the first is more relevant than with the second (stop-scharm), as 
only the former has all its constraints relaxed. Single top production has already been 
considered in the MRSSM via the process of tree level squark pair production, with the 
squarks decaying to top or lighter quarks and gauginos [17] . Our focus will be on single 
top production via top-up FCNC processes in SUSY QCD. 



2 Single top in SUSY-QCD 

2.1 Top quark FCNC couplings 

The additional quark flavour mixing possible in SUSY follows from a misalignment be- 
tween the squark and quark mass eigenstates. Taking the quark mass matrix to be 
diagonal, the mixing is encoded in the off-diagonal elements of the squark mass matrix. 
Contributions to the squark mass matrix come from the soft breaking mass terms as well 
as the superpotential, D-terms and soft breaking trilinear terms after electroweak symme- 
try breaking. Mixing between left and right (LR) labeled squarks is only possible through 
the electroweak breaking of the superpotential \x term and the soft trilinear terms. As 
both of these terms violate a U(l) R symmetry, the MRSSM cannot have LR squark 
mixing. Consequently, LL and RR squark mixing can be larger within the total mixing 
limits. To remain consistent with the analysis of Ref. [Hj we will consider only LL and 
RR mixing. The squark mixing matrix is then given by 



(M\) . . = (m 2 ).. + ml + M 2 cos(2/3) (T 3 - Q sin 2 9 W ) 



^3<$j3, (1' 



with q E {ul, d L , ur, d R }, m U3 = m t , m d3 = m h and i is the generation index of the quark 
mass eigenstate basis. The matrices fh 2 are soft SUSY breaking mass terms, which by 
SU(2) symmetry obey the relation fv? uL = fh 2 dL . The mass eigenstates of the squarks can 
be found by diagonalizing the above mass matrix: 



- q* (M 2 ) i3 q 3 = -q* a (u}Mp 9 ) J b , 



with 

?. = E( t/ S) *■ ( 3 ) 

i 

Here Uq is a unitary matrix and the index a denotes the squark mass eigenstate basis. 

The dominant top quark FCNC couplings in SUSY QCD are given by the one-loop 
bubble and triangle diagrams shown in Figure [Tj The flavour mixing matrices Uq enter 
through the standard quark-gluino-squark vertices that are also present in the MSSM. 
Specifically, via the Feynman rules: 

Qi- 9a- qla ■ - iK a V2g s (ul a )j A P a , 

Qi- 9a~ ft™ : _ iK, a V2gs(Uq a ) ai t A Pa, (4) 

with a G {L,R}, K a = {+1 : L, — 1 : R} and a = o{L -h- R). The amplitudes of these 
diagrams are vulnerable to a GIM-like suppression, as they each involve the sum over 
orthogonal elements of the unitary matrix U&: 

a£{L,R} a 

for i E {1,2}, where the function / represents the amplitudes dependence on the up 
squark masses. If the up squark masses are degenerate the amplitude will vanish. In 
general, the greater the mass difference (or splitting) in the squark sector, the larger this 
amplitude will be. Note that the UV divergent terms of these loop diagrams vanish by 
this mechanism. 
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Figure 1: The bubble and triangle diagrams contributing to the t-u-g FCNC. The arrows 
on the gluino lines are indicative of their Dirac fermion nature. 

2.2 Single top with Dirac gluinos 

As discussed in the introduction, the presence of Dirac as opposed to Majorana gluinos 
can suppress the SUSY-QCD box diagrams that contribute to neutral meson mixing 
phenomenology. The requirement for this suppression to occur is that the gluino, besides 
from being Dirac, is sufficiently heavier than the squarks in the loop. If this is the case, 
the strong constraints placed on first generation squark flavour mixing by meson mixing 
experiments are relaxed. In this paper we are concerned with single top production and 
are thus only interested in the mixing between the first and third squark generations. For 
simplicity, we parametrize the unitary matrices Uq by a single Euler angle 9 q for each 
q G {u L ,d L ,u R ,d R }. 

The production of a single top quark in the SM at tree level is accompanied by an 
additional lighter parton that eventually hadronizes into a jet. With the top-up FCNC 
discussed above, it is possible to produce exclusively single top quarks without this extra 
jet [12]. Our focus will nonetheless be on the production of a single top plus jet in SUSY- 
QCD, as it mimics the SM single top signal that is about to be extensively scrutinized at 
the upcoming LHC. Top-up FCNC processes that give single top plus jet are ug — > tg, 
gg — )■ tu, uq — >■ tq, uq — > tq and qq — > tu where q G {u, d, s, c, b}. In Figure |2] we present 
the diagrams that contribute to the process ug — > tg. 
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Figure 2: Feynman diagrams contributing to ug — > tg. The shaded t-u-g vertex represents 
the diagrams in Figure [TJ Cross-channel gluon diagrams are not shown. 

The relevant parameters for the SUSY-QCD FCNC couplings are the gluino mass m g , 
the squark mixing angles 9 u l, 9 uR and the 1st and 3rd generation squark masses m„^ a , 
iT^uRa- We make no attempt to calculate what regions of the parameter space satisfy 
the shifted mixing constraints beyond what was done by Ref. [Hj. Rather, we take a 
pragmatic approach by assigning a prototype point in parameter space and in turn vary 
the parameters with respect to it. We take this point to have a heavy gluino m g = 2 TeV, 
maximal squark mixing 9 u l = 9 u r = ir/4 and squark masses with a sizable mass splitting 
mj L = 400 GeV, m^ L = 1000 GeV. Right handed squark masses are set to 90% of their 
left handed counterparts. 



3 LHC Phenomenology 

3.1 Signal Features 

To calculate the partonic one loop squared amplitudes giving top plus jet we made use 
of the software packages FeynArts, FormCalc and Loop Tools [HI [191 EH EI]. The LHC 
hadronic cross sections were computed using the CTEQ6M PDF sets [22J , with the renor- 
malization and factorization scales set to the top quark mass and the LHC centre of mass 
energy taken to be 14 TeV. To stay within the jet detection limits at the LHC, a maximum 
pseudorapidity cut of 3 and a minimum transverse momentum cut of 20 GeV are placed 
on the outgoing parton accompanying the top quark. 

Of the single top plus jet processes possible via the top-up FCNC couplings, we found 
the process ug — > tg to be by far the most dominant at the LHC with a contribution 
of more than 80%. This can be attributed to the up quark density in the colliding 
protons, together with the exchange of t-channel gluons within its leading diagrams. At 
the prototype parameter point discussed in the previous section, ug — > tg has a hadronic 
cross section of 8.5 pb. This is in contrast to the runner-up process uu — > tu, which has a 
hadronic cross section of 404 fb. A model independent study of a t-u-g chromo-magnetic 
coupling reports similar results [23J. The dominance of ug —¥ tg could give an observable 
and characteristic asymmetry between single top and anti-top production at the LHC, as 
the process ug — > tg does not involve valence quarks. 

A key feature of the ug — > tg process is that the top quark is produced very forward, as 
shown in Figure 3(a) . This makes it distinct from SM single top, for which the differential 



cross section is peaked more centrally. For the semi-leptonic top decay this forwardness is 



carried over to the charged lepton, as shown in Figure 3(b) The detection of this signal 



at the LHC is therefore reliant on the forward efficiency of the detector. 
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Figure 3: The FCNC process ug —¥ tg, as well its semileptonic decay, have a distinct 
forward signature to their SM counterparts at the LHC. Figures (a) and (b) give the 
normalized pseudorapidity distributions for the top quark and charged lepton respectively. 
Only the parton accompanying the top has been cut on with 77 < 3 and pt > 20 GeV 



The one loop process ug — ¥ tg is dominated by the bubble and triangle t-u-g couplings 
shown in Figure [TJ A useful feature of these couplings is that their leading form factors 
have limited kinematic dependence and therefore vary little across the LHC phase space. 



This was tested numerically using FormCalc, with the triangle coupling taken to be on- 
shell to reduce the number of form factors from 16 to four. The triangle and bubble 
form factors that contribute most significantly are those involving the spinor structure 
I^Pl/r and Pl/r respectively It is thus possible to approximate the signal process by 
using constant tree level couplings with this spinor structure. 

3.2 Signal and background 

To assess the strength of our signal over the SM background at the LHC we used Mad- 
Graph/MadEvent to generate events [2H [25] . For the signal events an effective tree level 
t-u-g vertex with a constant form factor was used, as discussed in the previous section. 
Normalized against the full cross section as computed by FormCalc, the effective vertex 
was found to accurately reproduce the original kinematic distributions. Subsequently, the 
Madgraph DECAY package was used to decay the top plus jet signal to a charged lepton, 
neutrino, b-jet and jet. To simulate detector smearing effects, the four momenta of the 
final state quarks and gluons are scaled by a Gaussian distribution that is centred around 
their energy. The model parameters are again set to the prototype point discussed in 
Section E3 

The acceptance cuts for the signal and background are taken to be rj < 5, px > 15 GeV 
for positively charged leptons, rj < 3, px > 20 GeV for jets and px > 20 GeV for the missing 
transverse momentum. The large pseudorapidity cut for the lepton is in anticipation of 
its forwardness and will be addressed shortly. Valid events are required to have exactly 
two jets with AR > 3.5, of which one is b-tagged, and one positively charged lepton (with 
taus excluded). No restrictions are made on the number of negatively charged leptons. At 
the partonic level, i.e. without considering realistic showering or detector simulations, the 
main backgrounds to this final state at the LHC are SM single top, tt and W plus two jets. 
A b-tagging efficiency of 50% is assumed together with mistag rates of 0.5%, 1.5% and 
10% for the light quarks, gluons and charm quark respectively. To reduce the SM single 
top background the cut px < 75 GeV is placed on the non-b-tagged jet. This leaves W plus 
two jets as the dominant background in the forward region. To reduce this background, 
a cut is placed on the the top mass, 150 GeV < m t < 190 GeV, which is reconstructed 
from the four momenta of the charged lepton, b-tagged jet and neutrino. The unknown 
neutrino momentum is deduced from the total missing transverse momentum, where the 
constraint m\+ v = mw on the lepton-neutrino invariant mass is used to fix the longitudinal 
degree of freedom. We have also examined the effect of extra jets on these backgrounds 
and found their impact to be marginal under the selected cuts. The largest effect is from 
it plus one jet, which contributes an additional 10% to the non-dominant tt background. 

The lepton pseudorapidity distribution of our prototype signal and its SM background, 
after cuts, is shown in Figure |4j As expected, the signal is most promising with respect 
to the background in the forward regions of an LHC detector. Normalized against the 
central peak, there is a clear difference in shape between the presence and absence of the 
signal. A minimum pseudorapidity cut of 77 > 1.5 on the lepton gives a reasonable signal 
over background ratio. For instance, in the central bin the ratio is 2% whereas in the bin 
at rj ~ 2, where there is double the signal and almost half the background, the ratio is 8%. 
Of course, the upper bound of rj < 5 for the lepton is only to illustrate its forwardness. 
The LHC detectors ATLAS and CMS typically quote a pseudorapidity limit for muons 
between 2.5 and 3. With our selection criteria, and within the pseudorapidity window of 



1.5 < \r)\ < 2.7, the signal cross section is 100 fb and the ti(j), tj and Wjj background 
cross sections are 260, 320 and 600 fb respectively. An integrated luminosity of 10 fb -1 
therefore gives a statistical significance of S/yB = 9.6. 
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Figure 4: Pseudorapidity distribution of the final state lepton for the signal and SM 
background after cuts. As in in Figure 3(b) , the signal is seen to peak in the forward 
region. There is, therefore, a clear shape difference between the presence and absence of 
the signal when both are normalized against the central peak. 

Alternatively, our signal could be a candidate for the LHCb detector, which can mea- 
sure very forward leptons and jets. In this case, we would take 2 < rj < 5 for our 
acceptance cuts together with pt > 25 GeV for jets and pr > 15 GeV for leptons. Valid 
events are then defined as a charged lepton together with a b-tagged jet. The catch is 
that we cannot adequately reduce the tt background in the region where the signal peaks, 
as LHCb is blind to the second top quark. Beyond i] > 3 the signal over background and 
significance are again adequate. However, as the signal has already peaked, there is no 
longer a clear shape deviation. 

The results given so far have been at the prototype parameter point. We now briefly 
discuss how the relevant parameters each affect the signal cross section. A lighter gluino 
mass can significantly amplify the signal. Halving the gluino mass to m g = 1 TeV, for 
example, more than triples the cross section, whereas a heavier mass of m g = 3 TeV 
reduces it by 60%. We should keep in mind, however, that the mass of the gluino is 
related via experimental constraints to the level of flavour mixing allowed in the model. 
Decreasing the flavour mixing from its maximal value to 9ll — @rr — tt/8 has the effect 
of halving the cross section. The cross section scales with the squark masses in the same 
way as it does with the gluino mass. The size of the squark mass splitting, however, plays 
a crucial role for the cross section size due to the GIM-like mechanism present. Decreasing 



the mass splitting to 200 GeV reduces the cross section to a tenth of the prototype's cross 
section. 

We would like to reemphasize that the cross sections quoted in this paper were com- 
puted at leading order in QCD. A further quantitative analysis would involve a full next- 
to-leading order analysis, which could have a significant impact on our results. In partic- 
ular, due to the forwardness of the signal and the presence of soft p? cuts, the presence 
of initial state radiation could seriously affect the signal significance. 

4 Conclusions 

We have studied single top production at the LHC in a SUSY-QCD model with a heavy 
Dirac gluino. The presence of a heavy Dirac gluino allows for less constrained top-up 
FCNCs. In this scenario, the FCNC process ug — > tg was found to give the largest 
contribution. This can be attributed to the up quark density in the colliding protons 
and t-channel gluon exchange in the leading diagram. The key features of this signal 
are that the top quark is produced very forward and that it is asymmetric to its anti- 
top counterpart, as the latter lacks a valence quark. The semileptonic decay of the top 
in turns gives a forward charged lepton whose pseudorapidity distribution is distinct in 
shape from that of SM single top production. 

The signal was compared to its reducible SM background at the partonic level, without 
considering realistic showering or detector simulations. The prototype parameter point 
chosen was found to have a promising signal over background ratio and statistical signif- 
icance if a minimum pseudorapidity cut of 1.5 is placed on charged leptons. Moreover, a 
difference in shape between the presence and absence of the signal is clearly visible in the 
forward regions of the lepton pseudorapidity distribution. The feasibility of this signal is 
therefore dependent on the forward detection efficiency of muons at the LHC detectors. 
If a shape difference for positively charged forward leptons can be detected, it would be 
suggestive of the top-up FCNCs we have examined in this paper. 
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